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SUMMARY 

This paper presents a structural and aeroelastlc analysis of a large scale 
advanced turboprop rotor blade. This 8-blade rotor Is designed to operate at 
Mach 0.8 at an altitude of 35 000 ft. The blades are highly swept and twisted 
and of spar/shell construction. Due to the complexity of the blade geometry 
and Its high performance, It Is subjected to much higher loads and tends to be 
much less stable than conventional blades. Four specific analyses were con- 
ducted: (1) steady deflection; (2) natural frequencies and mode shapes; (3) 

steady stresses; and (4) aeroelastlc stability. State-of-the-art methods were 
used to analyze the blades Including a large deflection, finite element struc- 
tural analysis and an aeroelastlc analysis Including Interblade aerodynamic 
coupling (cascade effects). The study found the blade to be structurally sound 
and aeroelastlcal ly stable. However, It clearly Indicated that advanced turbo- 
prop blades are much less robust than conventional blades and must be analyzed 
and fabricated much more carefully In order to assure that they are structur- 
ally sound and aeroelastlcally stable. 


INTRODUCTION 

When propellers were replaced In commercial aircraft by turbojets, the 
overall propulsive efficiency was reduced, but the aircraft cruising speed was 
Increased. Since fuel was a small part of the overall operating cost, the 
Increased total passenger or cargo miles that could be provided by a given air- 
craft with the higher speeds more than made up for the Increased fuel consump- 
tion. Ever since the drastic rise In fuel costs during the 1970's, a major 
thrust within the aeronautics community has been to reduce aircraft fuel con- 
sumption. The development of low bypass ratio turbofan engines returned some 
of the lost efficiency without loss of performance and also resulted In more 
powerful engines. The subsequent development of high-bypass-ratio engines 
resulted In even further Improvements. As part of the overall effort to 
Increase propulsive efficiency even more, the NASA advanced turboprop program 
has brought the development of commercial aircraft propulsion full circle. 

The goal of this program Is to develop turboprops which can replace turbofans 
without any appreciable loss In aircraft performance, but with a significant 
gain In fuel economy. Analytical and experimental studies have shown that fuel 
savings of 15 to 30 percent are possible compared to today's turbofan engines 
(ref. 1). These turboprops have 8 to 10 low-aspect-ratio highly-swept blades 
which are twisted along the span and curved back about the axis of rotation. 

The forward speed of the aircraft Is about Mach 0.8 and the blade tip speed 
approaches 800 ft/sec. These turboprop blades will operate In a very compli- 
cated and severe aeromechanlcal environment. 


Designs such as these have been known for some time, but the means to 
analyze them and to build them did not exist. Today, with the development of 
high speed computers and sophisticated engineering software, advanced turbo- 
props can be analytically studied more thoroughly and more accurately. Also, 
with the development of advanced materials such as fiber reinforced composites, 
these advanced designs which meet structural and aerodynamic performance 
requirements are feasible and can be built. 

This paper presents a structural and aeroelastlc analysis of a prototype, 
large-scale, advanced turboprop blade which was designed by Hamilton Standard 
Division of United Technologies Corporation under contract to NASA. The blade 
(designated SR-7L and shown In fig. 1) Is part of a complete 8-blade 9-ft diam- 
eter rotor system scheduled to be built, tested and ultimately flown on a test 
bed aircraft. The results presented In this paper are from Independent anal- 
yses of SR-7L conducted at the NASA Lewis Research Center. 


Blade Description 

The SR-7L large scale advanced turboprop blade analyzed In this paper Is 
47 In long from the bottom of the shank to the tip. The airfoil section Is 
41 In long. The blade Is designed for a cruise Mach number of 0.8 at 35 000 ft 
altitude with a rotational blade tip speed of 800 ft/sec. The power loading 
under these conditions Is 32.0 shp/ft 2 . The blade Is of spar/shell construc- 
tion. The spar Is solid aluminum and blends Into an Integral aluminum shank. 
The spar Is 3.1 In thick at the base of the airfoil section and tapers to a 
thickness of less than 0.05 In near the blade tip. The shell Is of variable 
thickness and Is made from plies of woven fiberglass cloth. Each ply Is 
0.007 In thick. The blade twist and offset (relative to a radial line) also 
vary along the span and are shown In figure 2. 

An additional construction feature Is the nickel sheath along the leading 
edge of the blade. The purpose of the sheath Is to protect the blade against 
foreign object damage from small hard objects such as stones or hall. The 
thickness of the sheath varies from a maximum of 0.012 In near the blade lead- 
ing edge, to 0.006 In where It blends In with the shell. The overall complex- 
ity of the blade geometry and construction, along with the high power loading, 
tip speed and cruise Mach number, make the structural and aeroelastlc analyses 
much more difficult than for a conventional propeller blade. 


Analytical Procedure 

The objective of this study was to ascertain the structural stability and 
Integrity of the SR-7L blade design by analyzing four key areas of blade 
response: 


(1) Steady state displacements 

(2) Natural frequencies and mode shapes 

(3) Steady state stresses 

(4) Classical flutter 

In order to conduct these analyses, a variety of In-house and commercial 
computer codes were utilized. A brief description of the most significant 
codes follows: 
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(1) COBSTRAN - an In-house, special purpose composite mechanics program 
designed to produce composite-equivalent Input data sets for COSMIC NASTRAN and 
MSC/NASTRAN and to calculate Individual ply stresses from NASTRAN output data- 
sets. The composite mechanics In COBSTRAN Is fully generalized and contains 
general failure criteria. However, the code contains streamlined geometric and 
finite element model algorithms designed specifically for blade-like structures 
(ref. 2). 

(2) COSMIC NASTRAN - a large, general purpose finite element analysis 
computer program developed for NASA. 

(3) MSC/NASTRAN - a commercial version of NASTRAN with enhanced 
capabilities. 

(4) ASTROMIC - an In-house, fast-running classical flutter code designed 
for preliminary analysis and parametric studies of rotating blade- like struc- 
tures. The code uses subsonic and supersonic two-dimensional, cascade unsteady 
aerodynamics and models the structure as a straight, swept (or unswept) beam 
(ref. 3). 

(5) ACA/NASTRAN - a modification of the classical flutter analysis In 
COSMIC NASTRAN. This enhancement Incorporates the effects of cascade unsteady 
aerodynamics and blade sweep effects Into COSMIC NASTRAN which then uses a 
modal analysis to calculate aeroelastlc stability (ref. 4). 

(6) EZPLOT - an In-house graphics code streamlined for use with NASTRAN 
Input and output datasets. 

In addition to these codes, several special purpose pre- and post- 
processors were written to provide the proper Interfacing between the codes. 
Furthermore, several modifications had to be made to these major codes or to 
their solution procedures In order to fully accommodate the very complex 
geometry and structural behavior of SR-7L. 

The overall structural and aeroelastlc analyses followed the procedure 
shown In the flowchart In figure 3. Both analyses started with a geometric/ 
design database for SR-7L. A special preprocessor was used to extract detailed 
geometric data for the spar and shell In order to build a detailed composite- 
equivalent finite element model. While this model was being built, distributed 
beam properties derived from the database were used to build an ASTROMIC model 
In order to conduct preliminary aeroelastlc analyses. These analyses were used 
to estimate the aeroelastlc stability of SR-7L and guide the more refined aero- 
elastlc analyses conducted with ACA/NASTRAN. 

The structural and aeroelastlc analyses of SR-7L proceeded along essen- 
tially, Independent paths. The fundamental bridge between the two analyses 
was the detailed finite element model of the blade. 

The approach was to model the spar/shell construction by an equivalent 
laminated plate model. The shell, adhesive, spar, and shell filler material 
were taken to be separate layers and were combined to produce equivalent plate 
properties using the composite mechanics In COBSTRAN. This process, as It 
applies to SR-7L, will be described In greater detail In the next section. 

Other approaches, such as merging a separate three-dimensional spar model with 
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a hollow shell model or generating a combined three-dimensional spar/two- 
dlmenslonal shell model were considered. However, these two approaches greatly 
complicated the modeling process and would have necessitated major changes to 
the finite element modeling algorithms In COBSTRAN. Also, the approach adopted 
has proven very successful In modeling turbine engine fan blades Including 
superhybrid composite blades and blades with hollow cavities and composite 
Inlays (ref. 5). 

Once a COBSTRAN model was built, two NASTRAN Input datasets were gener- 
ated: one for MSC/NASTRAN, and one for COSMIC NASTRAN. Both versions of 

NASTRAN had to be utilized to analyze the blade structurally and aeroelas- 
tlcally due to the geometrically large deflections that can occur In thin, 
swept blades such as SR-7L. This has been observed experimentally (ref. 6) 
and predicted analytically for smaller aerodynamic models (ref. 7). While 
COSMIC NASTRAN cannot accurately account for this possibility, a modified 
application of the nonlinear displacement solution capability In MSC/NASTRAN 
(solution sequence No. 64) was used for detailed analysis of deflections, 
stresses and natural modes. 

On the other hand, a COSMIC NASTRAN Input dataset for the COBSTRAN model 
was needed for use In ACA/NASTRAN. Prior to conducting a detailed aeroelastlc 
analysis with this modified version of COSMIC NASTRAN, blade modal character- 
istics must be generated elsewhere In the code and supplied to the aeroelastlc 
analysis. COSMIC NASTRAN does not have the capability to rigorously account 
for large deflections. However, the code can Iterate to accurately determine 
the differential stiffness due to the softening or stiffening effects of non- 
follower type applied forces acting on the deflected blade. This approach Is 
not as accurate as the fully nonlinear analysis Incorporated In MSC/NASTRAN but 
can be acceptable for calculating reliable stability boundaries provided the 
discrepancies between the rigorous and approximate analyses are small. For all 
cases analyzed In this study, the agreement between the approximate and rigor- 
ous nonlinear analysis was sufficient to assure that good aeroelastlc analyses 
could be conducted using ACA/NASTRAN. This judgment was based on comparisons 
of the blade deflections and natural frequencies. 

Once the COBSTRAN model of SR-7L was built and the COSMIC NASTRAN struc- 
tural response was verified against the MSC/NASTRAN response, detailed struc- 
tural and aeroelastlc analyses proceeded Independently. 


STRUCTURAL ANALYSIS 
Finite Element Modeling 

The finite element mesh for SR-7L was developed from coordinate descrip- 
tions defining the blade spar and shell cross-sectional geometries at 27 
stations along the blade span. The midchord line of each cross section was 
determined from the aerodynamic profile data over a radial span of 10. 0 to 
54.0 In (as shown In fig. 4). Each midchord line was then represented by 26 
grid points defining coordinate values and thicknesses. 

COBSTRAN Is designed to accept 14 grid point definitions at each radial 
section along the blade span. Selection of grid points was made to establish a 
finer mesh along the leading and trailing edges of the spar (as shown In fig. 5). 
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This was done to minimize the effect of the elements representing a transition 
of properties between the all-shell nodes and the shell/spar nodes. This Is 
particularly necessary In C0BS1RAN as opposed to other programs, because of the 
special manner In which the program generates composite material properties at 
all grid points and then averages these properties among nodes defining an 
element In order to generate the material properties within each element. As 
such, the material properties of a grid point forming the boundary between the 
spar and the shell will contribute to the material properties of the elements 
on both sides of the boundary. 


Modeling Strategy 

The blade was designed as a fiber-glass/epoxy composite shell with an 
Inner spar of solid aluminum. The area within the shell and outboard of the 
spar was filled with foam material. Spar and shell were bonded with a layer of 
adhesive and the leading edge was covered with a sheath of nickel. A typical 
cross section Is shown In figure 6. 

The material properties at each node were determined by evaluating each 
node as a series of layers of fiber-glass/epoxy, adhesive, and foam or alumi- 
num with a symmetric stack-up through the thickness. The corresponding ani- 
sotropic stress-strain relations were calculated by laminate theory In which 
plane cross sections remain plane. Individual ply properties were calculated 
using micromechanics theory. Each node In the SR-7L model contained from 10 to 
66 ply layers with a different ply stack-up order depending on the location of 
the node on the blade. 

The COBSTRAN program Is designed so that ply layers may be defined to 
exist over only specified areas of the blade expressed by a percent span and 
percent chord range (as shown In fig. 7). Making use of this feature, and 
repeated Iterations through the COBSTRAN preprocessor, resulted In each of 326 
nodes being defined by their correct thickness and ply order. Ply order Is 
defined by sheath, shell, adhesive and foam or aluminum, as required, at each 
node location. 

The shank of the blade was modeled using a bar element. The transition 
between the triangular plate bending elements of the blade and the bar element 
representing the shank was modeled by the use of multipoint constraints. The 
root of the shank was fixed In 4 degrees of freedom, the three translations 
and rotation about the spanwlse axis. To represent the flexibility of the 
blade retention mechanism, rotation about both axes normal to span axis were 
represented by a rotational spring. 

A normal modes analysis (MSC/NASTRAN solution sequence No. 3) was first 
performed to determine the first four natural frequencies and mode shapes at 
zero rotational speed. 

Determination of steady state deflections, natural frequencies and mode 
shapes at selected rotational speeds was a more complex procedure. Because of 
the twist and sweep of the blade, steady state deflection calculations tend to 
be nonlinear. As such, the geometric nonlinear analysis capability available 
In MSC/NASTRAN was used (solution sequence No. 64). This procedure requires, 
as a minimum, two subcases to account for differential stiffness. Additional 
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subcases update the stiffness and mass matrices at each Iteration. The number 
of additional subcases Is problem dependent. For this analysis, eight addi- 
tional subcases were determined to be sufficient to establish convergence. 

A modification (DMAP alter) to the MSC/NASTRAN solution sequence No. 64 
was made to prevent numerical Instability. During the Iteration the stiffness 
matrix tended to become singular due to very small values of rotational stiff- 
ness about the z-axls. To prevent this, for the first seven Iterations the 
value of the stiffness matrix representing rotational stiffness about the z- 
axls was arbitrarily Increased by a comparatively small value of 0.001. This 
additional stiffness was then excluded from the eighth and final Iteration. 

When calculating the stiffness matrix change In a centrifugal force field, 
MSC/NASTRAN does not account for the softening effects resulting from displace- 
ments In the plane of rotation. To obtain the correct frequencies and mode 
shapes, this softening effect must be Included In the final stiffness matrix 
calculated In solution sequence No. 64. Therefore, a modification (DMAP alter) 
to the MSC/NASTRAN solution sequence was made to reduce the diagonal terms of 
the stiffness matrix associated with the translational degrees of freedom In 
the plane of rotation by a value of (m x fi 2 ). The final mass and stiffness 
matrices were stored by the MSC/NASTRAN database manager In temporary files for 
subsequent use by solution sequence No. 63, which Is the normal modes analysis 
with database files. 

The determination of ply stresses, strains and failure criteria Is made 
possible by the use of the postprocessing feature of the COBSTRAN program. 
COBSTRAN accepts the element stress output from MSC/NASTRAN and calculates the 
stresses and corresponding membrane forces and bending moments at each node In 
the structural coordinate system. From these loads, using laminate theory, the 
resultant stresses, strains and failure criteria are calculated for each ply 
layer at each node In the ply coordinate system. 


Steady State Deflections 

The steady state deflections of SR-7L at cruise are shown In figure 8. 

The blade Is loaded by centrifugal force and aerodynamic load, though the 
latter Is small compared to the former. Since the blade Is highly twisted, 
the contours are not shown In the global reference system. The contours 
represent deflections normal to a blade chord line at the three-fourths span 
position. As such, the deflections shown have very little edgewise component 
over the outer third of the blade where total deflections are the largest and 
are approximately normal to the blade In this region. The key feature of the 
deflection pattern Is that while the peak deflections are small compared to 
the overall blade dimensions, the deflections and gradients tend toward 
"moderate" values relative to the outer part of the blade. This tendency Is 
the result of the very thin blade construction, high centrifugal acceleration 
and built-in radial offset near the blade tip. Variations In any of these 
three factors can have a significant effect on the tip deflections. It Is 
possible for tip deflections to vary significantly from blade to blade as a 
result of deviations from design In the fabrication of the blade tip. This Is 
quite different from conventional straight propellers which are rugged along 
their entire length and essentially aligned with a radial vector so that 
centrifugal forces tend to restrict blade deflections, not Induce them. 
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Natural Frequencies and Mode Shapes 


The first four natural frequencies of SR-7L are shown on a Campbell dia- 
gram In figure 9. These span the frequencies of greatest concern with regard 
to flutter and forced response, although the first edgewise mode does not 
appear to contribute significantly to either. The Campbell diagram shows how 
the blade natural frequencies compare with engine order excitation frequencies 
as a function of the turboprop rotational speed. In order to avoid possibly 
severe dynamic response and subsequent fatigue damage, the critical blade 
natural frequencies should not be too close to an engine order excitation fre- 
quency (1-P, 2-P, etc.) near the turboprop operating speed. The "boxed-ln" 
regions about each engine order are the exclusion zones near the operating 
speed. The lower engine order excitations, In general, tend to cause stronger 
blade forces and, as such, their respective exclusion zones are larger than for 
higher engine orders. As can be seen, the SR-7L, as designed, meets all fre- 
quency placement requirements with the possible exception of the fourth engine 
order Interacting with the second bending mode. However, this Interaction 
typically would not be anticipated to be strong and probably would not cause 
this particular blade design to be rejected. 

Some special care must be taken In modeling blades like SR-7L for modal 
frequency analysis and In evaluating the results. Figure 10 shows the first 
four natural frequencies of SR-7L with and without a protective paint layer. 
This layer adds no stiffness to the blade but does add mass. Even though the 
layer Is only 0.01 In thick, It Increases the blade Inertial properties by up 
to 10 percent near the tip. From this figure, the first three frequencies are 
effected only slightly by the paint layer. However, the fourth mode, first 
torsion, Increased by about 3.5 percent when. the paint layer was excluded from 
the model. Since modal frequency exclusion zones are typically about 5 percent 
of the natural frequency for high modes, variations of a few percent can be 
significant. 

The mode shapes of the fully loaded SR-7L at the normal operating condi- 
tion are shown In figure 11. The contour lines are, as with steady displace- 
ments, shown for deflections normal to the chordllne at three-fourths span. 

The values of each contour line, subsequently, are properly scaled but rep- 
resent arbitrary units for each mode. 

The first mode appears as a "typical" blade first bending mode. The 
contour lines are essentially chordwlse and are fairly evenly spaced along 

most of the blade span. The second mode Is the first edgewise mode. Most of 

the motion In this mode Is In the chordwlse direction and, as such. Is not 
shown by the contour lines In this figure. However, since this mode Is not 
expected to contribute significantly to flutter or forced response, relative 
to the other modes shown, a detailed study of Its motion Is not needed for an 
analysis of SR-7L or similar blades. 

The third mode Is the second bending mode. The key feature of this mode 
Is that most of the blade motion Is In the outer third of the blade. Also, the 
contour lines tend to be In a chordwlse direction. In some cases, swept blades 
of this type will have modal contour lines that slant downward from the leading 

edge of the blade to the trailing edge. Such blades will then have bending- 

torsion coupling which tends to be aeroelastlcal ly destabilizing and, thus, 
lowers the flutter speed. For the SR-7L analyzed here, the absence of any sig- 
nificant bending-torsion coupling should enhance the calculated stability 
boundary. 
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The fourth mode Is a "typical" first torsion mode with a midchord mode 
line. The primary significance of this mode Is that It couples with the bend- 
ing mode aeroelastlcal ly to affect the stability of the blade. This form of 
coupling Is not characteristic of unswept blades, In general. 


Steady State Stresses 

The steady-state stresses In SR-7L at the cruise condition are shown for 
the aluminum spar and fiberglass shell In figures 12 and 13, respectively. The 
regions with the highest stress are shown In both figures. In the spar, two 
regions exist with stresses over 10 000 psl. One region Is on the pressure 
side of the blade near the bottom of the spar, and the other on the suction 
side of the blade toward the top of the spar. The lower region has a peak 
stress of 15 280 psl acting In a radial direction. The upper region has a peak 
stress of 15 050 psl, again In a radial direction. These stresses are domi- 
nantly bending, a result of the blade being offset from a radial line. 

The shell also has two regions of high stress In approximately the same 
relative locations as for the spar. The shell Is made from layers of cloth 
with a 0°, and -90° weave oriented In -30°, and *60° directions on the blade 
(as Indicated In fig. 13). The cloth was subsequently modeled as two uniaxial 
layers, one In the -30° direction and the other In the 60° direction. Both 
regions encompass stresses over 3000 psl, and In both regions the highest 
stresses are In the -30° ply. The peak stress In the upper region Is 6010 psl 
and the peak stress In the lower region Is 6194 psl. 

The peak stresses In both the spar and shell are relatively low and 
Indicate that the cyclic stresses are below the allowable vibratory stresses. 


C0NCLUSI0NS-STRUCTURAL ANALYSIS 

Based on the structural analysis conducted during this study, several 
general conclusions can be drawn regarding the SR-'/L advanced turboprop blade. 
In particular, and similar blades, In general: 

(1) Blade deflections at cruise are structurally acceptable. 

(2) Frequency margins relative to engine order excitations appear to be 
adequate (with the possible exception of the third mode, second bending, 
relative to the 4-P engine order). 

(3) Steady-state stresses at cruise are acceptable. 

(4) Blade deflections and critical vibration modes are dominated by 
motion of the comparatively thin outer portion of the blade (approximately the 
outer one-third). The results of the analysis and the actual behavior of the 
blade may be very sensitive to modeling and manufacturing detail In this 
portion of the blade. 
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AEROELASTIC ANALYSIS 


Aeroelastlc Methods 

As previously mentioned, two computer codes were used to determine the 
aeroelastlc characteristics of SR-7L. The major difference In these two codes 
Is the type of structural model. These are a beam model and a finite element 
model for ASTROMIC and ACA/NASTRAN, respectively. ASTROMIC does have the 
advantage that It Is easy and quick to use. ACA/NASTRAN Is a modification of 
COSMIC NASTRAN generated by Bell Aerospace under a NASA contract. The exten- 
sive capabilities of this code allow It to be used as either a design/analysis 
or a research tool. Because ACA/NASTRAN Is an extension of COSMIC NASTRAN, It 
can model complex structures. This capability makes the code well suited for 
modeling the structural configuration of a turboprop blade such as SR-7L. The 
disadvantage of using this code Is that It requires a large amount of computer 
time and memory, resulting In lengthy computer runs. Also, the ACA/NASTRAN 
user must be fairly knowledgable In finite element techniques to prepare the 
extensive Input data and to Interpret the results. 

The output from ACA/NASTRAN and ASTROMIC provide Information about both 
the dynamic and aeroelastlc characteristics of the turboprop blade. The 
dynamic characteristics supplied In the output are the natural frequencies and 
mode shapes. For the aeroelastlc part of the analysis the programs output the 
aerodynamic damping for each mode of vibration. The flutter boundary Is con- 
structed by finding the locations of zero aerodynamic damping as a function of 
axial Mach number and rotational speed. The following sections give brief 
descriptions of the two methods. 


Beam Method (ASTROMIC) 

The details of the derivation of the beam method for unswept blades are 
given In reference 3. This method was developed for pretwisted, nonuniform 
blades by using Hamilton's principle. The derivation of the equations has Its 
basis In the geometric nonlinear theory of elasticity In which elongations and 
shears are negligible compared to unity. A general expression for foreshorten- 
ing (axial shortening of the tension axis due to bending, torsion, and noncoin- 
cidence of the elastic and tension axes) Is explicitly used In the formulation. 
This method for unswept blades was modified In an approximate manner to account 
for blade sweep. For simplicity It was assumed that, although the blades are 
swept, the elastic axis Is straight. Only the component of centrifugal load 
along the blade axis Is considered. Also, for the flutter problem the blade 
Is assumed to be vibrating about Its undeformed position. The use of these 
assumptions can result In significant errors In the prediction of higher mode 
flutter characteristics. For example, It Is known that proper consideration of 
the steady state displacements can cause the pure torsion mode (nonrotating) to 
change Into a lower frequency, highly coupled bending-torsion mode. This modi- 
fied version of the reference 3 method also has the capability to consider 
blade mlstunlng, but this Is not considered herein. 

Both subsonic (ref. 8) and supersonic (ref. 9) two-dimensional unsteady 
cascade aerodynamic theories are used. The assumed relative flow Is the com- 
ponent normal to the elastic axes. The chord and stagger angle are also 
defined for sections normal to the elastic axis. The lift and moment are 
Integrated In a strlpwlse manner to give a quasl-three-dlmenslonal effect. 

These aerodynamic loads are corrected for sweep effects by using similarity 
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laws. This method, used in references 10 and 11 for an Isolated, nonrotating, 
swept wing. Involves modification of the two-dimensional lift and moment 
expressions for an unswept wing. The spanwlse component of flow Is neglected 
and similarities In the vertical velocity boundary conditions for swept and 
unswept wings are utilized. 

The space variable In the resulting coupled Integro-partlal differential 
equations of motion Is eliminated by using a modified Galerkln's method. The 
trial functions are the uniform beam mode shapes. For all results presented In 
this paper, two modes were used for each of three types of motion: bending In 
the plane of rotation, bending perpendicular to the plane of rotation, and tor- 
sion. The assumption of simple harmonic motion results In a generalized eigen- 
value problem which Is Iteratively solved to determine the flutter boundary. 
Further descriptions of this analytical method can be found In references 3 
and 12. 


Finite Element Method (ACA/NAS1RAN) 

The detailed description of this method Is given In references 4, 12 
and 13. For the reader's convenience, a brief description follows. The method 
consists of three steps: (1) a "differential stiffness" matrix Is determined 

from the steady state solution (by considering both centrifugal and steady 
aerodynamic loads and by using a Newton-Raphson Iteration to determine the 
equilibrium position); (2) by using the additional stiffness from step 1, a 
free vibration analysis Is performed to determine the modal characteristics; 
and (3) a modal flutter analysis Is performed. 

The two-dimensional, subsonic, cascade theory described In reference 14 Is 
used to calculate the unsteady aerodynamic loads. The theories of references 8 
and 14 are closely related and give similar results. A supersonic aerodynamic 
theory Including sweep effects has yet to be Incorporated Into this method. 

The subsonic unsteady aerodynamic loads are modified to account for sweep 
effects by considering only the component of flow normal to the local leading 
edge. Although this method can Include steady aerodynamic loads, they are not 
considered herein. 

In general, a blade Is a surface with multiple curvature. The method to 
form the aerodynamic forces for such a surface Is as follows. For each station 
a "mean surface" Is constructed between chord lines (which are normal to the 
leading edge and not parallel to each other) for this station and the adjacent 
spanwlse locations (Inboard and outboard). The modal translations normal to 
this surface are then used to construct the generalized aerodynamic force 
matrix. 

As In the beam method, a complex generalized eigenvalue problem results. 
The flutter boundary Is determined by repeated use of the "KE-Method" (ref. 15) 
for different axial Mach numbers and rotational speeds. Since this computer 
code presently lacks supersonic unsteady aerodynamics, the relative flow Mach 
number (perpendicular to the local leading edge) at all stations must be less 
than 1. This restricts the operating conditions at which this method can be 
used. 


10 



Previous Aeroelastlc Analyses 

Three other advanced turboprop designs (SR-5, SR-3 and SR-3C-X2) have been 
tested and analyzed for flutter. Since the conclusions and recommendations for 
SR- 7 L are related to the comparisons of these three experimental and anal /tlcal 
results, they will be briefly described. 

SR-5 and SR-3 are solid metal (T1) blades with different sweep distribu- 
tions. The SR-3 blade has slightly more sweep than SR-7L while the SR-5 blade 
has about 15° more sweep. SR-3 and SR-3C-X2 have the same geometric shape but 
consist of different materials, solid metal (T1) and graphite-epoxy, 
respectively. The ply orientations of SR-3C-X2 were specifically chosen to 
make the blade unstable within the capabilities of the experimental wind tunnel 
set up. The following sections will briefly describe the flutter analysis and 
experimental results for these three blades. 


SR-5 

A detailed description of the experimental results for SR 5 Is given In 
reference 12. Also Included In this reference are the results of a flutter 
analysis using the beam code which was described above. These results, 
assuming a constant sweep angle of 15°, are summarized In figure 14. As can 
be seen, there Is only fair agreement between the experimental and analytical 
results. The overall experimental flutter boundary agrees well with the cal- 
culated boundary for a setting angle of 61°. A comparison of flutter bounda- 
ries for different setting angles shows the same trend. That Is, Increasing 
setting angle Is a destabilizing effect. A comparison of flutter boundaries 
for similar setting angles shows only fair agreement between the experimental 
and calculated values. The calculated values always being somewhat conserva- 
tive. It should be pointed out that the beam method Is not expected to give 
precise values, since the structural model Is not complete. 

The SR-5 flutter analysis was repeated using the finite element method 
described previously. The finite element model for the blade Is shown In 
figure 15. The rotating frequencies and mode shapes were first calculated by 
performing a geometrically nonlinear steady state analysis and then using the 
final stiffness matrix In a subsequent linear free vibration analysis. The 
centrifugal softening terms were Included In both steps. The resulting 
Campbell diagram Is shown In figure 16. 

Also shown In the figure are experimentally determined natural frequen- 
cies (ref. 6). There Is good agreement between the predicted and measured 
frequencies for the first two modes. These frequencies and mode shapes were 
used In the modal flutter analysis, ACA/NASTRAN. The resulting flutter bound- 
ary Is shown In figure 17. As with the beam method, there Is good agreement 
between the overall experimental flutter boundary and the calculated boundary. 
However, when similar setting angles are compared, the agreement Is only fair. 
The beam and finite element method prediction are In good agreement and both 
are on the conservative side of the experimental boundary. 
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SR-3 


The SR-3 model blades were examined for flutter susceptibility both 
experimentally and analytically. Within the capabilities of the wind tunnel, 
11700 rpm and Mach 0.82, the blades were stable. The blade was analyzed using 
both the beam and finite element methods. The finite element model Is shown In 
figure 18. The beam analysis was conducted assuming a constant sweep angle of 
15°. In contrast to the experiment, both methods predicted flutter well within 
the test envelope. The predicted flutter boundaries are shown In figure 19. 

As can be seen, there Is fair agreement between these methods for low rota- 
tional speed and poor agreement for high speed. 

The beam method has the capability to consider both structural damping and 
mlstunlng (ref. 16). To determine the sensitivity of SR-3 to these two effects 
and possibly explain the discrepancy between analysis and experiment, addi- 
tional flutter boundaries were calculated considering mlstunlng and damping 
separately. The results are shown In figure 19. The use of 0.5 percent struc- 
tural damping resulted In a moderate effect on the flutter boundary, however, 
the blade was still predicted to flutter within the test envelope. The use of 
10 percent alternate mlstunlng of the bending frequencies resulted In a 
stronger effect. The effect Is strong enough to move the flutter boundary 
outside the test envelope. As a result of this study, It Is concluded that 
neglecting mlstunlng and/or structural damping may cause the flutter predic- 
tions to be overly conservative. 


SR-3C-X2 

In addition to the solid metal SR-3, there Is a family of blades with the 
same geometry but constructed of graphite/epoxy composite material. These 
blades are known as SR-3C . One set In this family, known as SR-3C-X2, was 
built-up with *22.5° ply orientations specifically designed to make the blade 
"flexible" and, hence, very susceptible to flutter. This set of blades was 
subjected to extensive wind tunnel testing and analysis. A subsequent publica- 
tion will describe the experiments and analyses In detail, however, they will 
be summarized herein. The finite element model used was the same as In 
figure 18 for SR-3 but the material properties were obtained using C0BSTRAN 
(ref. 2). 

The experimental and analytical results are shown In figure 20. As can 
be seen, there Is poor agreement. Once again, the analytical results are on 
the conservative side of the experimental flutter boundary. It should be 
mentioned that the analyses predicted the blade to be relatively Insensitive 
to setting angle. The reasons for the very conservative nature of the analysis 
Is an area of current study. Those areas being considered are mlstunlng, 
structural damping, three-dimensional aerodynamic effects (tip unloading and 
effective sweep), and transonic aerodynamic effects. 


SR-7L FLUT1ER ANALYSIS 
ASTR0MIC Analysis 

Since the SR-7L blade has variable sweep and has spar/shell composite 
structure, equivalent beam properties must be estimated. Therefore, the 
C0BSTRAN finite element model was used to calibrate the ASTR0MIC SR-7L Input 
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data. Adjustment was then made to ASTROMIC material properties In order to 
obtain good agreement In the first bending, second bending, and the first tor- 
sion modes (fig. 21). This approach was only good for these lower modes and 
starts to break down In the higher modes. ASTROMIC does a poor job of predict- 
ing the edgewise mode frequencies. This Is seen In the large discrepancy 
between the codes for the first edgewise mode shown. The disagreement Is 
attributed to the fact that ASTROMIC does not Include shear deformation and 
rotary Inertia. 

Based on the geometric finite element model, the assumed straight elastic 
axis of SR-7L was swept 10°. This constant sweep angle resulted In the effec- 
tive Mach number distribution shown In figure 22. Processing of the SR-7L 
Input data results In the ASTROMIC modal damping predictions shown In 
figure 23. This figure displays the relationship between damping and flight 
Mach number at a specific operating condition (1700 rpm, sea level). The curve 
for this mode shows a sharp Instability at 0.67 Mach number which would not be 
easily eliminated with Increased damping. ASTROMIC predicted Instability In 
the first mode only. 

To assess the effect of altitude on the stability of SR-7L, the same 
ASTROMIC analysis was repeated at a constant rotational speed (1700 rpm) and 
various altitudes. This produced the ASTROMIC flutter boundary represented by 
the dashed line In figure 24. Increasing the altitude appears to have a stabi- 
lizing effect on SR-7L as can be seen by the movement of the flutter boundary 
to the right as altitude becomes higher. Ideally, no flutter prediction should 
occur to the left of the flight profile, but as can be seen In the figure, the 
flutter boundary crosses the flight envelope at 0.80 Mach number. This was not 
cause for concern In that ASTROMIC has been shown to underestimate the flutter 
speeds of SR-5. The reasons for this conservatism will be addressed later In 
this discussion. 


ACA/NASTRAN Analysis 

The effective Mach number used by ACA/NASTRAN differs considerably from 
that used In ASTROMIC. The difference comes from the way that the aerodynamic 
loads are corrected for the blade sweep. While ASTROMIC Is limited to a 
constant sweep angle (straight elastic axis), ACA/NASTRAN can account for the 
actual variation In the blade sweep. The effective Mach number In ACA/NASTRAN 
Is taken as perpendicular to the leading edge and will thus change not only 
with the rotational velocity component but also with the sweep angle. The 
effective Mach number distribution for SR-7L for ACA/NASTRAN, along with that 
of ASTROMIC Is shown In figure 22. 

ACA/NASTRAN Is restricted In that It does not have transonic or supersonic 
capability and, therefore, cannot handle relative Mach numbers of 0.95 or 
greater. With this limitation and the specified rotational speed of 1700 rpm, 
the maximum flight Mach number Is fixed at 0.88. Figure 25 presents damping as 
a function of flight Mach number for three modes at four flight speeds (0.40, 
0.60, 0.80 and 0.88 Mach numbers) and at a constant rotational speed of 
1700 rpm (sea level). Neither mode one or three approaches Instability at 0.88 
Mach number. Mode two Is a predominantly edgewise bending mode for which the 
aerodynamic loads are always small. 
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Referring once again to the flight envelope depicted In figure 24, the 
ACA/NASTRAN aeroelastlc analysis of SR-7L was performed at 1700 rpm with 
various altitude conditions. No Instability could be found up through the 
maximum Mach number (M = 0.88) capability of the code. 

The results generated by both ASTR0M1C and ACA/NASTRAN for SR-7L are 
believed to be conservative. First, damping was not Included In the analysis. 
Disregarding damping Is conservative since both material and mechanical (such 
as root friction) damping are known to have a stabilizing effect on flutter. 
Second, blade mlstunlng was neglected. Results from previous studies on the 
effect of mlstunlng (ref. 17) Indicate that a moderate amount of mlstunlng can 
alleviate flutter problems. Third, both ASTROMIC and ACA/NASTRAN use two- 
dimensional aerodynamic theories which assume that there are an Infinite number 
of blade cross sections out the span. This assumption does not allow for the 
Inclusion of tip unloading effects which would be Included were a full three- 
dimensional aerodynamic theory used. These three dimensional effects have been 
shown to move the flutter boundary to a higher operating condition. Neglecting 
the three-dimensional effects results In lower operating speed flutter predic- 
tions. The final reason for conservatism Is the limitations of the computer 
codes In their abilities to realistically simulate the transonic effects. Both 
codes use a less than adequate estimation of the transonic effects which 
appears to result In conservative flutter predictions. 


CONCLUSION-AEROELASTIC ANALYSIS 

1. As previously described, the agreement between analytical and experi- 
mental flutter boundaries has ranged from fair to poor. However, the analyt- 
ical predictions have always been on the conservative side; that Is, the 
predicted flutter Mach number has always been less than that measured In the 
wind tunnel. The suspected reasons for this conservatism have been discussed. 

2. SR-7L Is predicted to be stable over the entire flight envelope by the 
finite element method. The more approximate beam method predicts this blade to 
be stable over most of the flight envelope. A small region In the upper right 
hand corner Is unstable. In light of these predictions and the knowledge of 
the conservative nature of the analytical predictions, the SR-7L Is judged to 
be a low risk aeroelastlcal ly . 


SUMMARY OF RESULTS 

Four separate analyses were performed on the SR-7L advanced turboprop 
blade: (1) a steady state deflection analysis; (2) a modal analysis; (3) a 

steady state stress analysis; and (4) an aeroelastlc stability analysis. The 
first three analyses, forming a complete structural analysis, were based on a 
large deflection finite element analysis. The aeroelastlc analysis used two 
Independent approaches: (1) an approximate beam analysis, and (2) a refined 

finite element analysis. This study found the SR-7L blade to be structurally 
sound and aeroelastlcal ly stable. However, there are generic characteristics 
of blades of this type which can cause less carefully designed blades to be 
structurally and aeroelastlcally unacceptable. 

First, most of steady state deflections and vibratory motion occur In the 
outer third of the blade which Is very thin. This part of the blade Is also 
heavily loaded and, as such, the blade tip undergoes geometrically moderate 

14 



deflections. Also, due to sweep the blade can have a degree of vibratory 
pitch-flap coupling which is detrimental to stability. 

Because the blade is so thin, small deviations from design specifications 
during fabrication can have a significant effect on the overall blade charac- 
teristics. This was demonstrated by comparing the blade natural frequencies 
with and without a thin protective paint layer. This layer caused one of the 
critical natural frequencies to change by 3.6 percent. Considering that reso- 
nance margins are only about 5 percent, a variation of this size can be signi- 
ficant. As a result, special attention must be given to designing, analyzing 
and fabricating the outer portion of blades like SR-7L. 

It was also found during this study that the flutter boundaries predicted 
by the approximate and refined analyses differed significantly, though the 
overall analysis indicated that the blade would be stable over its entire 
proposed flight envelope. While the approximate beam analysis predicted some 
flutter at high Mach numbers and altitudes, the more refined analysis, with a 
better structural and aerodynamic model, predicted no flutter. Subsequently, 
the approximate analysis proved to be a good tool for qualitatively examining 
the aeroelastic stability of SR-7L but is not adequate for a detailed study of 
such a complex blade. For that matter, even the refined analysis has defi- 
ciencies in that three-dimensional aerodynamic tip effects, transonic effects 
and structural damping are not Included. However, these effects tend to 
Improve stability, and as such, the aeroelastic analysis conducted can be 
considered conservative. 

In conclusion, while this study found the SR-7L advanced turboprop blade 
to be structurally sound and aeroelastical ly stable, special care must be taken 
when designing, analyzing and fabricating blades of this type. Also, though 
the analysis methods used in this study proved to be adequate for SR-7L, there 
is room for considerable Improvement in methodology for structurally and aero- 
elastically analyzing advanced turboprop blades, as well as in construction 
and fabrication concepts. 
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Figure 3. - Analytical procedure. 
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Figure 4 - Cross-sectional profile. Blade geometry provided as separate 
profiles for spar outline and shell mid-thickness. 
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Figure 6. - Typical blade cross-section. 



Figure 7. - COBSTRAN laminated plate model. 
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Figure 8. - Steady deflections at cruise conditions 
(1700 rpm with air loads). 
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(c) Third mode, m.6 Hz, (d) Fourth mode, 16a 9 Hz, 

second bending. first torsion. 

Figure 11. - Natural mode shapes. (Contour lines are arbitrary units normal 
to the blade chord at the 3/4 span position, } 
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Figure 12. - Regions of highest spar stress. 




Figure 14 - SR-5 flutter data. 
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Figure 19. - SR-3 flutter data. 
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Figure 20. - SR-3C-X2 flutter boundary. 




Figure 21. - Comparison of ASTROMIC and ACA/NASTRAN 
natural frequencies. 
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Figure 22. - Effective Mach number distortion. 
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